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ABSTRACT 


A  method  for  resolving  the  mesoscale  distribution  of  precipitation  over  mountainous 
regions  is  developed  using  synoptic-scale  ancdyses  and  topograph}'  with  a  40  km 
resolution.  The  method  consists  of  calculating  the  saturated  moishue  flux  at  a  level 
assmned  to  be  the  cloud  base  and  correlating  this  quantity  with  observed  precipitation. 
Geireral  correlation  coefficients  ranged  from  -0.2  under  conditions  poorly  represented  by 
the  svnoptic  surface  analysis  to  -0.55  when  these  analyses  are  more  representative  of  actual 
forcing.  When  applied  regionaUv,  correlation  improved  dramatically,  suggesting  potential 
incorporation  into  precipitation  forecasts  over  mountainous  areas.  Over  wider  regions  of 
less  defined  topography,  correlation  is  weaker  but  is  stiU  an  improvenrent  over  the  current 
synoptic-scale  precipitation  forecast  products. 
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viii 


I.  introduction 


The  enhancement  of  synoptic-scale  precipitation  by  topography  is  well  understood, 
and  has  been  routinely  verified.  Fine-mesh  mesoscale  models  have  proven  particularly 
accurate  in  reflecting  orographic  influences  on  large-scale  precipitation  events  (Colton 
1976,  Collier  1976).  Simpler  estimates  of  this  enhancement  factor  have  been  shown  to 
be  almost  exclusively  functions  of  regional  topography  combined  w'ith  low-level  wind 
fields,  synoptic-scale  vertical  motion,  and  a  knowledge  of  the  low-level  relative  humidity 
(Browning  et  aJ.  1975.  HUl  et  al.  I'^l^l ). 

Thus  far.  only  detailed  mesoscale  models  have  been  able  to  reproduce  the  small- 
scale  orographic  enhancement  process.  With  synoptic-scale  precipitation  products 
showing  little  or  no  regional  structure,  and  mesoscale  forcing  fields  generally  unavailable, 
forecasters  rely  mostly  on  local  knowledge  and  vague  statistical  multipliers  to  predict 
areas  of  heavy  rainfall.  Consequently,  the  ability  to  quantify  local  orographic  effects  and 
then  associated  distribution  of  precipitation  based  solely  on  synoptically  derived  wind  and 
moisture  fields  would  be  of  particular  use  to  the  forecaster. 

In  this  study,  the  above  mentioned  quantities  of  subsynoptic-scale  topography  and 
synoptic-scale  wind  and  moisture  data  are  combined  in  an  effort  to  reproduce  areas  of 
enhanced  precipitation  during  a  .series  of  eastern  Pacific  extratropical  cyclones  which 
passed  over  California  during  the  period  1200  IJTC  27  looi  1200  UTC  04 

Ntn.rCll  I  I  ritl,  t-..’  .-,1  »  I  ,  ,  x  p  |x\  ||if> 
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National  Meteorological  Center’s  (NMC)  final  analyses,  produces  a  distribution  of  the 
total  vertical  moisture  flux.  This  vertical  moisture  flux  estimate  is  correlated  vvith  hourly 
precipitation  data  collected  by  the  National  Oceanic  and  Atmospheric  Administration 
(NOAA),  covering  244  reporting  stations  throughout  California  for  the  same  period,  to 
produce  an  equivalent  rainfall  amount. 

The  overall  intent  of  this  study  is  to  aid  forecasters  in  identifying  areas  most  likely 
to  experience  orographic  enhancement  of  precipitation  for  a  given  synoptic  event  and 
provide  a  quantitative  precipitation  estimate.  This  procedure  may  be  of  use  in  estimating 
seasonal  snowpack  levels  in  otherwise  inaccessible  regions,  as  weU  as  providing  near  term 
flood  potential  for  specific  drainage  basins,  while  only  requiring  synoptic-scale 
information.  More  specific  objectives  of  this  study  are  to: 

-  Assess  the  degree  of  accuracy  to  which  precipitation  distribution  can  be  represented 
using  mesoscale  topography  and  large-scale  wind  and  moisture  fields. 

*  Develop  a  simple  quantity,  calculated  locally  from  model  output  products,  to 
represent  the  enhancement  of  precipitation  due  to  topography. 

*  Present  specific  synoptic-scale  conditions  and/or  regional  topographic  characteristics 
which  either  promote  or  discourage  the  use  of  this  parameter. 

-  Determine  the  level  of  improvement,  if  any,  the  employment  of  this  parameter 
shows  over  current  precipitation  forecast  products. 

Following  a  discussion  of  background  work  done  in  the  field  of  orographic 
precipitation  enhancement,  an  outline  of  the  dynamical  basis  for  the  vertical  moisture  flux 

a,'!  a  precipitation  pritrM'n'^n  i  <.  J'l--  •■•in  I  >  ■  . I  tlv 

inpnf  field"  Next,  fhi':  pm  '  .-,,1 1 .1  .it.-. I  p..  •  .p,t  M.'n  oV^t.iin  ;t 
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statistical  fit  assigning  specific  precipitation  accumulations  to  the  calculated  values 
followed  by  an  examination  of  the  factors  which  affect  correlation.  Results  obtained  by 
applying  the  precipitation  parameter  both  generally  and  regionally  are  examined,  including 
a  comparison  to  current  precipitation  forecast  aids  (specifically,  the  Nested  Grid  Model 
(NGM)  12  hour  precipitation  forecast).  Conclusions  focus  on  the  usefulness  of  this 


technique  as  a  forecast  aid. 


IJ  BACKGROUND 


It  is  certainly  no  secret  that  mountainous  areas  often  receive  greater  amounts  of 
precipitation  than  do  upwind  or  downwind  low-lying  regions.  However,  it  was  not  until 
Bergeron  (1965)  introduced  his  "seeder  cloud  -  feeder  cloud"  concept  that  a  specific 
mechanism  for  topographically  enhanced  precipitation  became  accepted.  As  described  by 
Raddatz  and  Khandekar  (1979)  this  process  consists  of  low-level  (feeder)  clouds  forming 
in  upslope  areas  which  "wash  out"  when  seeded  by  raindrops  falling  from  synoptic -scale 
(seeder)  clouds  above.  The  idea  is  that  the  feeder  cloud  is  constantly  replenished  and  will 
continue  to  enhance  precipitation  in  proportion  to  its  depth,  provided  seeding  from  above 
is  uninterrupted.  Storebo  (1976)  and  Bader  and  Roach  (1977)  were  able  to  confirm 
Bergeron's  theories  by  results  which  related  a  10%  increase  in  precipitation  to  a  1000  m 
feeder  cloud  depth.  Interestingly.  Browning  etal.  (1974)  suggested  that  in  some  instances 
(namely,  cases  of  mid-level  potential  instability),  orography  alone  could  be  responsible 
for  generating  precipitation,  even  if  synoptic-scale  seeding  was  not  available.  According 
to  Browning  et  al.  (1975),  condensation  over  mountainous  regions  is  roughly  proportional 
to  the  low-level  wind  velocity  component  normal  to  the  topographic  barrier  combined 
with  synoptic-scale  vertical  motion  and  moisture  quantities.  This  relationship  was 
generally  verified  by  the  mesoscale  numerical  model  developed  by  Collier  (1975). 

As  mentioned  in  the  introduction.  Cohon  (1^76)  achieve^  Temnrkablv  accurate 
results  using  a  fine-ni‘:'';ii  (  i  '  l.-m  rTi<i)  n-i-xi.  I  , . '■  I  .....I  I..  -  I  applied 
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to  the  Sierra  Nevada  of  California.  A  key  assumption  made  by  Colton  was  to  define  the 
orographic  precipitation  as  a  quantity  dependent  on  the  presence  of  synoptic-scale 
precipitation.  As  such,  his  scheme  is  onJy  able  to  enhance  synoptic-scale  precipitation, 
and  cannot  force  precipitation  based  solely  on  orography. 

A  later  study  by  Hill  er  al.  (1981 )  further  confirmed  the  seeder-feeder  enhancement 
theor>  .  They  additionally  showed  that  if  the  total  precipitation  rate  exceeded  a  given 
thieshold  value  (0.5  mm  hr  '),  the  topographic  enhancement  factor  was  not  dependent  on 
the  background  precipitation  rate  In  short,  as  long  as  the  background  rainfall  rate  was 
sufficientlv  high,  enhancement  was  simply  proportional  to  the  low-level  moist  flow 
interacting  with  topography,  and  vvas  not  a  function  of  the  amount  of  seeding  from  above. 
This  carries  the  important  implication  that  orographically  induced  precipitation  can  simply 
be  added  to  synoptic-scale  rainfall  to  obtain  a  total,  rather  than  the  orographic  term  being 
a  function  of  the  synoptic-scale  tenn.  This  relationship  provides  the  basis  for  deriving 
the  precipitation  parameter  described  in  the  following  chapter. 

More  recent  studies  have  attempted  to  modify  feeder-seeder  theory  to  account  for 
additional  effects  including  stratification  and  wind  drift.  While  showing  the  effect  of 
stratification  to  be  negligible  as  it  relates  to  total  washout,  Carruthers  and  Choularton 
(1983)  did  show  wind  drift  effects,  which  were  not  accounted  for  in  models  such  as  Bader 
and  Roach  (1977).  to  be  detrimental  to  both  local  maxima  and  total  enhancement  by 
"advecting"  orographic  precipitation  to  downwind  locations.  Additionally,  this  study 

backed  up  ^he  ifle^  i'’-.  iiiii  .  ,  :•  ■  '  ■  '  ■  ••  ..i- . .  ■  ..vi'  prn  I'  lu  <>( 

seeding  t.ate  (a.Iso  co'''tt:tdi' tr'i\  t..  tI-..’  :itv'  i,,,.  i,  (  i’’  .  ,,,,,,  i,  i. 
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Most  recently.  Stow,  Bradley  and  Gray  (1991)  conducted  a  study  with  results 
largely  consistent  with  seeder-feeder  theories,  and  emphasized  the  significant  differences 
between  the  previous  enhancement  studies.  They  regarded  the  work  so  far  conducted  as 
inconsistent  in  terms  of  aling  with  the  effects  of  the  low-level  wind  speed,  seeding  rate, 
and  barrier  height.  Moreover,  they  felt  that  although  the  basic  physics  of  topographic 
enhancement  are  well  understood,  refinement  is  required  to  handle  the  above 
inconsistencies.  The  present  study  is  not  concerned  with  any  of  these  "refinements",  but 
rather  with  using  the  existitig  theoiy  to  derive  a  simplified  method  of  representing 
mesoscale-resolution  topographically  enhanced  precipitation  from  exclusively  synoptic- 
scale  data. 

It  should  be  noted  that  topography  plays  a  larger  role  in  weather  than  simply 
enhancement  of  precipitation.  In  fact  it  has  been  shown  that  topography  can  influence 
the  low-  to  mid-level  flow  in  such  a  manner  as  to  completely  reposition  large-scale 
pressure  centers  (Sawyer  1959),  as  well  as  greatly  altering  the  small-scale  wind  patterns 
(Danard  1977).  An  attempt  to  combine  small-scale  deformation  of  the  wind  field  due  to 
topography  with  established  enhancement  theory  was  undertaken  by  Anderson  and  Nilsson 
(1990)  and  had  only  moderate  success,  mostly  owing  to  large  errors  in  obtaining 
precipitation  measurements  in  drifting  snow.  Any  deformation  of  the  wind  field  by 
topography  will  be  ignored  in  this  thesis,  in  terms  of  calculating  our  precipitation 
paiameter.  As  such,  a  main  objective  of  this  study  is  to  determine  the  feasibility  of 

niesosca^*'  -U  . '.  i,  ,  ;*  tlv/  .'v  tn;d 

mesoscalc  wind  firh!~ 
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in  PROCEDURE 


A.  THEORY 

Preceding  studies  have  estimated  orographic  enhancement  of  precipitation  by 
determining  the  vertical  extent  of  a  feeder  cloud  which  forms  in  upslope  regions  based 
on  low-level  winds  and  wet  bulb  potential  temperature.  This  cloud  depth  can  be  directly 
related  to  orographically  induced  precipitation  provided  seeding  from  above  (synoptic- 
scale  rain)  is  available.  Assuming  the  "washout"  efficiency  is  properly  handled,  numerical 
models  can  easily  detennine  the  enlrancement  factor  based  on  calculated  feeder  cloud 
depth  (Raddatz  and  Khandekar  1979,  Bader  and  Roach  1977). 

Since  the  purpose  of  this  study  is  to  provide  an  easily-calculable  means  of 
estimating  topographical  enhancement  of  precipitation,  a  slightly  different  approach  than 
a  numerical  model  simulation  is  taken.  This  estimate  of  topographic  precpitation 
enhancement  is  based  on  a  "ground-up"  quantity  reflecting  the  total  saturated  moisture 
flux  at  a  given  height  (assumed  to  be  the  cloud  base),  (wq)^.  Conceptual  restraints  are 
imposed  such  that  any  upward  saturated  moisture  flux  at  this  level  is  condensed,  and  not 
allowed  to  either  be  redistributed  into  higher  levels  of  the  atmosphere,  or  result  in 
supersaturation.  Essentially,  any  such  flux  is  completely  "washed-out"  as  precipitation. 
Since  this  flux  is  representative  of  precipitation  only  under  saturated  conditions,  a 

threshold  700  mb  hnmnlit-.  ■  •.  I  . i-  tf> 

represent  saturati'^n  a'  tli*-  nopti.  t  >1'  n.  :,p,.!,.  .i  |..  th.  iv,  ''-'Id 
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B.  DERIVATION  OF  PRECIPITATION  PARAMETER 


where  DqJDr  represents  rate  of  change  of  liquid  water  content  following  the  air  parcel 

trajectory’.  C  is  used  to  represent  the  condensation  rate  while  A  is  the  rate  of  washout 
from  abr've.  If  follows  that  DqJDt  is  sitnply  the  total  precipitation  minus  evaporation 


wliere  F  and  E  stand  for  precipitation  and  evaporation,  respectively.  Following  Holton 
(1979,  pp  40-44), 

dq^  dq^ 
dr  dr 


^9.  ^<1,  ^<1. 

u - +v - +(0 - 

dx  dy  dp 


represents  the  total  time  derivative  of  q^,  in  terms  of  the  local  tendency  and  advective 

components.  Rearrangement  in  vector  notation  and  use  of  the  isobaric  form  of  the 
continuity  equation  with  (2),  gives 
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of 


(3) 


Changes  in  the  saturation  specific  humidity  {q  )  result  from  temperature  and  pressure 
changes,  which  are  largest  due  to  vertical  advection.  Thus,  the  local  tendency  dqjdf  is 


smaller  than  the  flux  divergence  V  qU  and  will  be  ignored  in  this  calculation.  The  P-E 
term  is  the  net  condensation,  which  is  equivalent  to  the  rainfall  when  integrated  over  the 
depth  of  the  atmosphere  Tlie  net  condensation  is  written  simply  P  in  the  following 
equations.  We  are  left  with 


-P=ViqU). 


(4) 


as  a  first-order  approximation  to  the  net  precipitation.  Separating  (4)  into  horizontal  and 
vertical  components  yields 


-P 


diqu)  ^  d{qv)  ^  9(^(0) 

d.\  dy  dp 


This  can  be  integrated  from  the  surface  to  the  top  of  the  atmosphere  H,  to  give 
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d{qu)  ^  d{qv) 


(5) 


Ph 

where  (^(OL  represents  total  moisture  flux  at  H,  which  is  zero.  The  term  - 

'  M 

P« 


represents  the  rainfall.  R.  and  the  right-hand  side  of  equation  (5)  gives  the  contribution 

Ph 

from  both  topographic  uplift.  ,  and  synoptic-scale  forcing,  * ^i.qU)]dp .  The 

Po 


synoptic -scale  forcing  is  simply  the  net  horizontal  flux  divergence  of  saturated  air  into  the 
colvtmn.  Because  q^  is  largest  at  low  levels,  the  net  moisture  flux  is  dominated  by  the 

lower  part  of  the  atmosphere.  Consequently,  this  term  is  approximated  by  the  integrated 
moisture  flux  convergence  from  the  surface  to  700  mb.  The  approximate  rainfall  is  then 
given  by 

700 

p, 

with  R  representing  net  precipitati'^n  in  the  eohimn  tprin-;  npi"' ’-’nr  on  th*'  rieht- 

hand  side  of  the  equat’o”  --.-i-il-,  olMph^- •!  P'-'.lni.i  i- ,.,.i  moi^me  data 
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combined  with  surface  topography.  The  40  km  resolution  topography  used  in  this  study 


is  shown  in  Figure  I.  Specificallv,  the  term  (^to)  is  given  as  q  [« — — ], 

dx  dy 

with  u,  V  and  h  representing  horizontal  wind  velocity  components  and  surface  slope, 
respectively,  while  refers  to  the  surface  moisture  value. 

The  surface  forcing  iqw)  .  and  ultimately  (^to)  .  was  calculated  using  the  1000  mb 

wind  field  and  topograpbv  in  meters,  with  the  1000  mb  winds  treated  as  crudely 
representative  of  synoptic-scale  surface  winds  for  this  purpose.  The  resultant  surface- 
forced  moisture  flus  is  in  units  of  .y/A.y  rm/j.  An  approximate  conversion  to  pressure 
coordinates  is  made  by  applying  the  relation  q\\  --p^(X>q . 

7t*f» 

The  synoptic  scale  tenn  j\^{qU)]dp).  is  calculated  using  a  centered  finite 

r„ 


difference  scheme  to  compute  the  moisture  flux  divergence  at  levels  specified  by  the  wind 
aitd  moisture  field  inputs.  Summations  over  the  desired  vertical  distance  yield  the  net 
moisture  flux  divergence  in  the  surface  to  700  mb  layer,  which  is  the  vertical  moisture 
flux  qco  through  the  700  mb  level.  This  is  the  quantity  mO/.  referred  to  in  the  Results 

chapter.  Figures  2  and  3  show  sample  contoured  surface  vertical  velocity  and  <7(0  flux 

fields  obtained  fr<^m  tti«'-;f'  ,  Tti-  l•■I . .  >i  ' . .  ’•  r-  'n  r  ii-nf-  ' 

a  »  iin. .  p  .#,.!  .......ij,,,..  [■'  '  ■  if'  f'n'.ii rp  tins 
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quantity  is  of  the  same  order  of  magnitude  as  the  synopticaJly  derived  forcing  in  Figure 
3.  FinaJly.  the  total  moisture  flux  at  700  mb,  (hereafter  referred  to  as  wqf),  is 

obtained  by  simply  adding  the  topographic  and  synoptic-scale  components,  utilizing  the 
rough  equivadency  of  pbars  '  and  cms  '  as  described  above  (Figure  4).  Thus  both  wO/and 

vi  qfare  in  units  of  gkg  1  pbars  I  .  This  unit  definition  applies  to  all  figures  showing  vif>/ 
or  fields  unless  otherwise  noted. 

Note  that  this  quantity  is  not  in  units  of  rainfall.  The  actual  associated 

rainfall  is  obtained  by  6rst  employing  a  15%  700  mb  relative  humidity  threshold  to  thefo)^).,,,,^ 

field  to  ensure  saturation  This  field  is  then  correlated  with  observed  12  hour 
accumulated  precipitation  centered  on  the  analysis  time.  A  linear  regression  is  calculated, 
yielding  a  conversion  factor  based  on  the  slope  and  y-intercept,  allowing (a>7),j^ 

quantities  to  represent  12  hour  precipitation  totals.  For  the  overall  correlation,  no 
objective  limitations  are  imposed  other  than  the  relative  humidity  threshold.  Regional 
constraints  are  imposed,  however,  when  examining  the  effectiveness  of  this  technique  in 
more  localized  mountainous  regions. 

Perhaps  the  most  obvious  limitation  of  using  this  quantity  to  estimate  precipitation 
is  its  exclusion  of  sub-grid  scale  (i.e  .  convective)  precipitation  processes.  While  no  effort 
is  made  to  quantify  the  effect  of  topography  on  such  precipitation,  these  conditions  are 

identified  using  satellite  (macerw  nnH  n.-t.-.i  v  M-i'-.-p-.ni-  •' . i- -i-  niv) 

section. 
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Figure  1.  Topography  (40  km  resolution,  in  meters) 


C.  DATA 

The  saturated  moisture  Oua  quantity  used  to  represent  precipitation  is  calculated 
from  NMC  final  Global  Data  Assimilation  System  (GDAS)  analysis  fields,  which  were 
then  statisticallv  fitted  h';’nr1v  pr^ripitpitirtn  i.-.i  '  ■  r to  \  a  Tlio 
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- i - ^ - il.  ^  '.i 

Figure  2.  Smoothed  orographically  forced  surface  vertical  velocity  (cm/s).  1200 
UTC  01  March  1991.  Negative  values  correspond  with  upward  velocities. 

NMC  GDAS  analyses  include  twelve  pressure  levels  from  1000  to  50  mb  (as  described 

by  Kanamitsu  1989).  For  this  study,  the  2.5'’  NMC  GDAS  grid  is  interpolated  to  the  40 

km  grid  on  which  the  regional  topography  was  derived.  The  wind  and  moisture  fields 

at  the  in0('.  8,^^  an^  ■  v- r.  >  ...  u.  ,i.  .|  n-druMt'' 


flu.x.  \\qf. 


Figure  ?.  SynopticoMy  forced  700  mb  moisture  (wO/0.  1200  LITC  01  March 
19^1  Negative  valiies  correspond  with  upward  flux. 
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Figwrf  4.  Total  7(K)  mb  moisture  flux  {wqf).  1200  UTC  01  March  1991. 
Negative  values  corresponci  with  upward  flux. 
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Specific  rairtfall  totals  obtained  through  hoiiriy  precipitation  data  publications 
(NOAA  for  the  period  '^f  study  were  used  to  correlate  precipitation  with  the  ^qf 

values  NOAA's  precipitation  totals  consist  of  hourly  rain  gauge  observations  made  at 
244  locations  within  California.  For  our  purposes,  the  observations  have  been  grouped 
in  six-hour  increments  exactly  preceding  and  following  the  NMC  analysis  time.  This 
gives  12  hour  precipitation  totals  centered  on  the  NMC  anadysis  time  which  are  then 
cf'rrelated  witli  the  corresponding  uy/f  values. 
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IV.  SYNOPTIC  ANALYSIS 


A.  OVERVIEW 

The  series  of  synoptic-scale  events  that  occurred  during  the  period  of  this  study 
were  produced  by  the  development  of  a  500  mb  trough  off  the  Northern  California  Coast. 
The  system  center  migrated  northeastward,  over  the  period,  and  supported  a  series  of 
short-wave  troughs  and  associated  surface  cyclones  (Figures  5  and  6). 

An  important  factor  in  the  synoptic  evolution  of  this  period  is  the  movement  of  the 
primary  jet  stream.  Approximately  midway  into  the  period  (1200  UTC  02  March  1991), 
the.  jet  repositioned  itself  over  Central  California,  from  the  upper  Baja  Peninsula  where 
it  was  located  early  in  the  period.  This  shift  is  evident  in  the  300  mb  isotach  analyses 
representative  of  the  jet  position  at  1200  UTC  28  February  1991  and  0000  UTC  03  March 
1991  (Figure  7).  Thie  jet  axis  shifts  from  a  westerly  to  a  more  southwesterly  orientation 
coincident  with  its  latitudinal  shift.  As  a  result,  the  weather  affecting  California  was 
dominated  by  highly  baroclinic,  fast  moving  cyclones  in  zonal  flow  with  the  surface 
center  well  poleward  of  the  main  jet  early  in  the  period,  and  weaker,  less  baroclinic 
systems  late  in  the  period,  with  only  a  brief  transitional  time  in  between. 

As  a  consequence  of  the  poleward  migration  of  the  jet,  incoming  weather  systems 
affected  more  Northern  sections  of  California  late  in  the  period,  as  evidenced  by  the  500 

mb  height  rjjes  ob."e'‘'’‘''i  '^•  <^1  '  :ti,i . ;•  .  ■  . i  ,  m,'-  ^  --f  ilv- 
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rain  is  concentrated  in  Southern  California  early  in  the  period,  with  virtually  no  rain 
observed  in  Southern  California  by  the  end  of  the  period. 
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Figure  5.  500  mb  heights  (meters,  left  column)  and  surface  pressure,  including 
frontal  positions  (mb,  right  column)  for  (from  top)  1200  UTC  27  February  1991, 
1200  UTC  28  February  1991  and  1200  LITC  01  March  1991 


FigHr€  6.  500  mb  heights  (meters,  left  column)  and  surface  pressure,  including 
frontal  positions  (mb.  right  column)  for  (from  top)  1200  UTC  02  March  1991,  1200 
irrc  03  March  1991  and  1200  UTC  04  March 
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B.  CYCLONE  CHARACTERTSTICS 


1.  0000  UTC  29  February  1991  -  1200  UTC  02  March  1991 

As  mentioned,  cyclones  affecting  the  weather  in  California  during  this  period 
developed  on  the  cold,  cyclonic  side  of  the  main  jet,  as  seen  in  the  GOES  infrared 
imagery  at  1200  UTC  28  February  1991  (Figure  8)  over  the  Eastern  Pacific  and  the 
corresponding  mean  sea  level  pressure  and  1000-500  mb  thickness  analyses  (Figure  9). 
This  particulai  image  and  analysis  are  indicative  of  the  general  nature  of  the  cyclones 
during  this  period.  The  strong  surface  thickness  gradient  along  30° N  as  well  as  the  low- 
level  instability  upstream  of  the  main  frontal  cloud  bands  on  the  GOES  imagery  indicate 
the  strength  of  the  baroclinic  systems  passing  through  California.  These  systems  had 
significant  pre-frontal  cloud  and  precipitation  bands  associated  with  them.  In  general,  this 
ty]?e  of  pattern  is  similar  to  that  in  which  Collier  (1975)  applied  his  mesoscale 
enliancement  model  over  North  Wales,  thus  providing  a  good  basis  for  comparison. 

Cyclones  of  this  nature  are  favorable  for  observing  orographically  enhanced 
precipitation,  largely  due  to  the  forcing  being  well  represented  by  the  large-scale  analyses 
and  the  precipitation  primarily  associated  with  the  occluded  front.  However,  post-frontal 
convection  does  pose  a  problem.  While  Collier  (1975)  applied  his  technique  to  rain 
associated  with  both  cold  and  wann  frontal  passages,  Bader  and  Roach  (1977)  tried  to 
focus  only  on  the  mostly  stable  precipitation  produced  in  the  warm  sector  of  cyclones. 
Browning  er  af.  (1974)  earlier  concluded  that  oroEmphicaHv  induc  t  rrecipitafion  w-as 
most  evident  in  the  pmi  >■.  f,,i  1,.^  i.  .  i .  .  :  i  and 
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after  cold  frontal  passage  StiJI.  tbei.r  results  indicated  orographic  effects  could  be 
quantified  even  during  and  after  the  passage  of  the  cold  front.  Results  obtained  by  Collier 
(197?)  showed  less  degradation  of  quantifiable  orographic  enhancement  during  cold 
frontal  passage,  as  compared  to  warm  sector  precipitation,  than  did  the  results  presented 
by  Browning  et  al.  (1974).  Collier  (1975)  placed  more  emphasis  on  distinguishing 
between  fast  moving  baroclinic  systems  in  which  mesoscale  forcing  is  minimized,  and 
slower,  less  baroclinic  features  as  far  as  the  observability  of  orographic  precipitation 
enhancement.  Since  the  cvcioties  io  this  period  of  the  present  study  can  be  characterized 
as  occluded  surface  features  with  mostly  stable  precipitation  prior  to  frontal  passage,  and 
some  post-fiontal  convection,  orographic  effects  should  be  well  represented.  This  study 
does  not  examine  in  detail  the  effect  of  mesoscale  post-frontal  convection  on  the  overall 
correlation  between  uv//  and  observed  precipitation. 
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Figure  7.  300  mb  height?  (solid,  meters)  and  isotachs  (dashed,  knots)  for  1200  UTC 
28  Febmary  1991  (top)  and  1200  UTC  02  March  1991  (bottom). 
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2.  1200  IJTC  02  M?rch  1/>0I  -  1200  UTC  04  March  1991 


The  generaJ  nature  of  the  three  weather  systems  affecting  California  during  this 
period  is  represented  in  Figures  10  and  I  i  The  thickness  gradient  has  shifted  northward 
off  the  California  coast  and  the  amplitude  of  the  pre-f;untal  thickness  ridge  is  less  than 
in  the  earlier  cyclones.  In  addition  there  is  a  lack  of  post-frontal  instability  depicted  on 

the  infrared  image.  An  exception  to  this  pattern  is  the  area  of  cloudiness  at  approximately  35° 

N.  150'’  W.  As  with  the  systems  preceding  it,  this  particular  system  becomes  cut-off 
from  the  cold  air  within  which  h  developed,  and  evolves  mostly  in  response  to  the  strong 
southwesterly  flow  aloft  This  system  (which  reaches  California  at  approximately  1200 
Lire  04  March  1991).  like  the  two  seen  ahead  of  it.  has  little  representation  of  a 
thickness  ridge  ahead  of  the  front  as  it  reaches  land.  Without  the  pre-frontal  baroclinic 
forcing,  the  rainfall  associated  with  this  type  of  system  is  less  dominated  by  stable  pre¬ 
frontal  precipitation  than  cyclones  in  the  earlier  period.  Due  to  the  uninterrupted  supply 
of  warm,  moist  low-  and  mid-level  flow  from  the  southwest  during  this  period,  incoming 
weather  systems  produced  less  organized  and  more  convective-type  precipitation.  This 
type  of  pattern  was  specifically  avoided  in  terms  of  quantifying  orographically  induced 
precipitation  by  Collier  (1975).  since  the  forcing  responsible  for  the  distribution  of 
precipitation  was  difficult  to  represent,  even  on  the  mesoscale.  Consequently,  this  period 
of  the  present  study  may  not  be  adequately  represented  by  the  synoptic-scale  and 
orographic  effects  on 
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3-  Summarv 


The  above  synoptic  discussion  suggests  two  cases  that  produced  the  observed 
rainfall  distributions; 

-  Case  (i)  -  strongly  baroclmic  cyclones  with  a  well  defined  pre-frontal  thickness 
ridge. 

*  Case  (iii  -  subtropically  influenced  convective  complexes  with  little  or  no  pre- 
f'ontaJ  thickness  ridge 

The  first  case  tends  to  dominate  the  winter  rainy  seaion  over  California, 
although,  the  second  case  is  not  uncommon.  Consequently,  results  from  the  two  periods 
in  this  study  cover  contrasting  synoptic  conditions  which  produce  rain  in  California. 
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V.  RESULTS 


A.  OVERALL 
L  Statistics 

Table  1  shows  the  correlation,  slope,  and  y-intercept  from  a  linear  regression 
analysis  of  and  observed  12  hour  precipitation  for  four  representative  analysis  times. 
The  correlation  values  run  from  a  virtually  non-correlative 

-0.19  to  a  marginally  respectable  -0.55,  as  shown  by  the  representative  scatter  plots  in 
Figures  12  and  13.  The  correlation  values  are  grouped  at  -0.2  and  -0.5  for  the  later  and 
earlier  synoptic  periods,  respectively.  This  suggests  that  the  synoptic  and  mesoscale 
factors  discussed  in  the  previous  chapter  predetermine  the  success  of  this  technique  to 
be  either  moderately  correlative  (-0.5)  or  essentially  non-correlative  (-0.2),  with  little  in 
between.  Averaging  the  slope  and  y-intercept  of  all  "correlative"  times  (Table  1  shows 
just  representative  times)  yields  -0.010  and  +35  respectively.  Therefore, 

(-(m-^x0.01]+35=T 

equates  calculated  w  q/  values  with  expected  12  hour  accumulated  precipitation  (P),  in 
inches.  Note  that  all  values  of  P  than  zero  represent  a  precipitation  value  of  zero. 

As  will  be  shown  in  a  later  section,  correlation  is  greatly  unproved  in 
mountainous  subregions  as  compared  m  tb^-  v  •  a';--  technique  is 

intended  to  be  "nnive''ar  n’  it  i.-  .•>  ••  i  <  ;ile  paratnciets 
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and  topography,  (6)  will  be  applied  in  all  sub-regions,  rather  than  deriving  specific 
statistics  for  individual  sub-regions. 

TABLE  1  LINEAR  REGRESSION  STATISTICS  {WQF  VS.  12  HOUR 
PRECIPITATION) 


CORRELATION 

COEEHCIENT 

SLOPE 

y- 

INTERCEPT 

(INCHES) 

0000  UTC  28 
FEB  1991 

-.55 

-.013 

+0.50 

1200  Lire  01 

MARCH  1901 

-.54 

-.010 

+0.35 

0000  Lire  03 
MARCH  1991 

-.003 

+0.40 

1200  UTC  04 
MARCH  1091 

-.10 

-.003 

+1.00 

One  interesting  statistical  result  is  the  fairly  consistent  y-intercept  value  for  the 
linear  regression  of  about  +3S  (average  of  six  "correlative"  times).  This  indicates  that  the 
precipitation  threshold  (the  .v-intercept)  generally  corresponds  with  shghtly  downward 
saturated  moisture  flux  values.  Though  the  smoothing  process  tends  to  distribute  upward 
flux  values  into  lee  slope  regions,  it  is  still  likely  out-weighed  by  the  "spillover"  eflfect 
observed  by  Colton  (1975).  In  this.  spUlover  process,  upward  vertical  velocities  are 
advected  over  the  summit  region,  effectively  "mixing"  upward  vertical  velocity  air  into 
lee  regions  This  process  is  not  accounted  for  in  the  Miscalculation,  which  sees  all  lee 
regions  as  having  d':>\^n^■^ mrl  in  ti>-  o.  . . •  .  .i 
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FiguT‘€  9.  Surface  pressure  (mb,  solid)  and  1000-500  mb  thickness  (meters, 
dashed).  1200  UTC  28  Februm  1991. 

Additionally,  horizontal  wind  drift,  as  described  by  Collier  (1975),  may  be  producing 
artificially  high  precipitation 

observations  across  summit  regions,  especially  under  conditions  of  strong  low-level  flow 
with  a  substantial  peipendicular-to-barrier  component. 

The  other  consistent  statistical  feature  is  the  slope  value  of  either  about  -0.01 
for  the  correlated  times,  or  about  -().003  for  times  showing  less  correlation.  This  suggests 
a  "flatter"  precipitation  distribution  during  non-correlative  times,  with  non-orographic 

piocesse"^  motp  Hc'mitvmi  m  o.,.  p,.  ■  .p.i  .>•  — 
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2.  ^  <lf  Pij:fribiifi<>n 

Although  the  linear  regression  statistics  computed  above  have  a  relatively  low 
correlation,  comparison  of  horizontal  distributions  of  the  inferred  precipitation  and 
observed  precipitation  reveals  qualitatively  correct  distribution.  In  fact,  this  comparison 
reveal^  a  significant  improvement  over  NGM  model  forecast  precipitation. 

a.  0000  rjTC  28  February  1991  -  1200  UTC  02  March  1991 

Figures  14  and  15  depict  the  statewide  precipitation  distributions  produced 
bv  it'O/  (the  synoptiC'Scale  forcing  described  in  Chapter  ID).  NMC's  Nested  Grid  Model 
(NGM^  forecast  precipitation  for  1200  LITC  28  Febniary  1991,  acnial  observations  and 
u  q/'  Of  note  here  is  the  similarity  between  the  NGM  forecast  precipitation  and  the  wO/ 
field,  panicularly  the  maximum  over  Southern  California.  This  implies  that  the  synoptic 
forcing  we  are  using  in  our  calculations  of  ^^^qf  bears  rough  equivalency  to  the  NGM 
precip'^ation.  This  time  is  especially  informative  since  the  12  hour  700  mb  vertical 
velocity  forecast  verifies  particularly  well  compared  to  the  kinematically  calculated 
analvzed  700  mb  vertical  velocity  (not  shown).  It  should  be  mentioned  that  vrO/ showed 
virtually  no  correlation  with  observed  precipitation  throughout  the  study. 

In  general,  we  see  a  significantly  greater  amount  of  detail  within  the 
precipitation  field  produced  by  ^cqf.  compared  to  either  M)f  or  the  NGM  forecast.  Also 
note  the  "false”  maxima  produced  in  the  observed  field.  This  is  mainly  a  result  of  the 
multiquadrant  contouring  teclmique  e.xtrapolating  ba-d<;  function';  into  roj-ipns  of  sparse 

da(a  ( i  e . .  r  ^oiaco  t  1 1  :in  I  ,M  •  ill--  . >■  •  !••  •  -  ' ■  •  •  i  '.  n .  ■  ■  i m ■  w/jf 
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Figure  10.  GOES  infrared  imagery  1200  UTC  02  March  1991. 


distribution  reflects  the  observations  reasonably  well,  particularly  the  maxima  observed 
over  the  mountains  of  Southern  California  and  the  Sierra  Nevada.  Note  that  correlations 
are  performed  between  the  gridded  wqf  fields  and  station  locations  (rather  than  gridded 
precipitation).  It  should  be  noted  that  overall  correlation  is  generally  constant  from  0000 
Lire  28  February  199)  to  0000  UTC  01  March  1991,  with  the  exception  of  time  0000 
UTC  01  March  1991,  which  ha.";  a  correlation  coefficient  of  -0.21.  This  drop  in 
correlation  is  most  likely  in  response  to  the  low-level  convection  seen  in  Figure  8  located 
between  the  main  i-nn.!  :  T  i,;-;  i-mm.  ••  ■  ',:i  aiulj.- , 
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post-frontal  precipitation  distribution  is  determined  less  by  orography  than  convective 
processes. 

b.  1200  UTC  02  March  1991  -  1200  04  March  1991 

Figures  16  and  17  depict  the  same  parameters  as  Figures  14  and  15 
except  for  1200  LfTC  04  March  1991.  which  is  representative  of  the  second  pteriod.  Here 
vt  <7/is  generally  a  less  accurate  representation  of  observed  precipitation,  particularly  along 
the  Central  California  Coast.  StUl.  orographic  maxima  produced  by  wqf  over  the 
mountains  of  Northern  California  and  the  Northern  Sierra  Nevada  are  reasonably 
representative  of  observed  maxuna.  Though  the  corresponding  NGM  precipitation 
forecast  show’s  good  areal  coverage,  it  is  devoid  of  any  of  the  above  structure.  In  this 
case  the  NGM  vertical  velocity  forecast,  on  which  precipitation  is  largely  based,  shows 
a  strong  area  of  upward  motion  located  over  Northern  California,  whereas  the  verifying 
kinematically  calculated  analysis  shows  this  area  of  upward  motion  to  remain  off-shore 
(not  shown).  As  such,  the  NGM  precipitation  field  owes  its  accuracy  more  to  chance  than 
to  improved  dynamical  representation.  It  seems  that  neither  the  NGM  forecast  or  the 
GDAS  analysis  fields  are  able  to  pick  up  on  the  forcing  responsible  for  the  locally  heavy 
rainfall  in  tlris  period.  Referring  back  to  Figure  6,  we  see  an  analyzed  "front"  off  the 
Northern  California  coast  at  this  time,  with  most  of  the  state  within  the  warm  sector.  A 
random  check  of  station  observations  supports  the  strong  southerly  flow  depicted  by  the 
GDAS  analyses.  Convective  lifting  would,  in  this  ca,':e.  be  the  best  explanation  for  the 
lack  of  upward  verf'’;;’'  •  u-,  ii>.  .  <  •••  . .  i.  \  ! 
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Figure  11.  Surface  pressure  (mb,  solid)  and  1000-500  mb  thickness  (meters, 
dashed).  1200  UTC  02  March  1991. 


earlier,  this  period  is  dominated  by  only  marginally  baroclinic  features  in  the  moist 
southwesterly  flow  aloft.  As  such,  convection  is  likely  to  be  equally  important  in  both 
the  pre-  and  post-frontal  periods.  Though  convective  processes  are  parameterized  in  the 
NGM  precipitation  forecasts,  they  are  not  accounted  for  at  all  in  the  calculation  of  wqf. 
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3.  Correlation  L'niiting  Fsietor? 
o.  Synoptic  Conditions 

As  described  by  Collier  (1975),  synoptic  conditions  favorable  for 
observing  orographic  enhancement  of  precipitation  are  characterized  by  fast  moving 
baroclinic  systems.  The  thinking  is  that  slower  moving,  or  non-baroclinic  systems  allow 
more  non-orographic  and/or  mesoscale  mechanisms  to  influence  rainfall  distribution.  This 
is  not  to  say  that  topography  does  not  influence  mesoscale  structure.  Rather,  since  only 
the  large-scale  structure  is  considered,  topography  is  only  allowed  to  influence  the 
synoptic -scale  fields.  As  mentioned  previously,  the  period  prior  to  1200  UTC  02  March 
1991  was  characterized  by  baroclinic  features  well  represented  by  the  synoptic  fields, 
whereas  non-orographic  and/or  a  mesoscale  processes  dominated  after  this  time.  The 
statistical  correlations  support  this  characterization  of  the  dominant  precipitation  type 
during  the  periods.  However,  as  seen  in  the  horizontal  distribution  of  wqf  during  the 
second  period,  the  qualitative  distribution  of  precipitation  is  realistic  in  the  second  period. 
This  suggests  that  the  amount  of  precipitation  is  poorly  represented  in  more  convective- 
type  periods,  but  the  general  distribution  is  still  strongly  influenced  by  the  topography. 

b.  Smoothing 

Aside  from  subgrid-scale  effects,  smoothing  of  the  topographic  m'  used 
to  calculate  seems  to  be  most  detrimental  to  correlation  values.  As  previously 

discu.",':ed.  "meso.'^cal'r  nphi'-  t<>t  •  < . ; . (  0(U) 

mb  winds  over  40  Pni  Te'-'-lntiott  -.  •( .  '•  p  ’  . to  i.o 
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Figure  1,2.  Observed  precipitation  (vertical  axis,  inches)  vs.  v^’qf  (horizontal  axis) 
and  calculated  linear  regression  for  0000  UTC  28  February  1991  (top)  and  1200 
UTC  01  March  1991  (bottotn). 
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compatible  with  the  synoptic-scah'  forcing  As  a  result,  stations  upwind  of  topographic 
barriers,  such  as  Fresno  (station  FKES,  Figure  18)  are  generally  assigned  values  far 
in  excess  of  the  actual  calculated  forcing.  Note  (Figure  18)  that  the  smoothed  value 
for  FRES  (-15)  is  nearly  tw'ice  that  of  the  unsmoothed  value  (-8). 
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Location 


Certain  locations,  such  as  coastal  areas,  are  more  likely  to  have  higher 
actual  relative  humidities  than  inland  areas,  which  are  analyzed  sis  equal  on  the  synoptic- 
scale.  As  a  result,  more  of  the  vertical  moisture  flux  \^‘qf  is  likely  to  be  converted  to 
precipitation  in  coastal  areas  than  inland  areas,  where  a  larger  amount  of  moisture  flux 
(wqf)  may  be  required  to  maintain  or  reach  saturation.  This  relationship  is  evident  in 
comparing  the  linear  regression  corresponding  to  the  San  Bemaidino  -  San  Gabriel 
subregion  (subregions  are  defined  in  the  following  sections)  to  that  of  the  Southern  Sieaa 
Nevada  (Figure  19).  Note  that  a  wqfvzlw  of  -100  represents  approximately  a  2.75"  12 
hour  precipitation  accumulation  according  to  the  linear  regression  derived  for  the  San 
Bemardino/San  Gabriel  subregion,  but  less  than  1.0"  based  on  the  Southern  Sierra  Nevada 
subregion's  statistics. 

<f.  Other  Considerations 

Obviously,  numerous  other  considerations  affect  the  correlation  of 
precipitation  and  700  mb  saturated  moisture  flux.  Redistribution  of  moisture  above  700 
mb.  quantification  of  evaporation,  and  wind  channeling  induced  by  the  topography  are  a 
few  examples.  Though  important,  such  processes  caiuiot  be  easily  accounted  for  on  the 
synoptic  scale.  In  fact,  these  factors  will  require  more  specific  consideration  if  wqf 
cannot  adequately  represent  the  forcing  patterns  which  lead  to  precipitation  distributions 
in  mountainous  regions. 
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Figure  13.  Observed  precipitation  (vertical  axis,  inches)  vs.  wqf  (horizontal  axis) 
and  calculated  linear  regression  for  0000  UTC  03  March  1991  (top)  and  1200  UTC 
04  March  1991  (bottom) 


Figure  14.  NGM  12  hour  precipitation  forecast,  valid  1200  UTC  28  February 
19^1  (mm,  top)  and  (bottom)  for  1200  LITC  28  February  1991. 


Figure  15.  ^^  qf  for  1200  LnrC  28  February  1991  and  observed  precipitation 
(inches,  bottom),  0600  UTC  28  Febniaiy  1991  -  1800  28  Febniary  1991. 
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Figure  1,6.  NGM  12  hour  precipitation  forecast  valid  1200  UTC  04  March  91 
(mm.  top)  a,nd  mO/ (bottom)  for  12(X)  UTC  04  March  1991. 
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Figure  17.  Wqf  (top)  for  12UU  U  IL  04  March  IVVl  and  observed  precipitation 
(inches,  bottom),  0600  LH'C  04  March  1991  -  1800  UTC  04  March  1991. 


Figure  1?.  Unsmoothcd  (top)  and  smoothed  (bottom)  orographic  ally  forced 
vertical  velocity  (cm/s).  QQW  LrTC  01  March  1991. 
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R.  RFGTONAI. 

1.  Approach 

''■Vhlle  the  general  results  are  indicative  of  the  overall  value  of  the  wqf 
parameter  in  representing  precipitation  over  a  wide  region,  it  is  perhaps  more  important 
to  determine  the  extent  to  which  the  correlation  can  be  improved  when  applied  to  a 
specific  topographic  region.  In  this  way,  wqf  may  be  able  to  identify  specific  regions  or 
even  drainage  basins  likely  to  receive  heavy  precipitation  under  certain  forecast 
conditions.  Foi  this  reason,  the  specific  perfonnance  of  wqf  in  regions  of  w'ell  defined 
mpographv  is  examined  to  detenu ine  the  ability  of  wqf  to  represent  precipitation  in 
specific  regions.  In  effect,  the  estent  to  which  synoptically  derived  low-level  winds 
interacting  with  topographs  can  represent  the  mesoscale  distribution  of  precipitation 
within  mountainous  regions  w  ill  be  detennined.  Little  or  no  improvement  of  the  regional 
correlation  coefficient  over  the  general  case  would  suggest  limited  usefubess  of  the 
synoptic-scale  winds  in  determining  orographic  forcing.  An  increase  in  correlation  would, 
however,  indicate  that  a  significant  component  of  the  precipitation  forcing  is  due  to  the 
large-scale  vvind  field  interacting  with  regional  topography,  and  that  this  quantity  is  more 
easily  isolated  in  regions  of  well  defined  topography.  Specific  correlation  coefficients  are 
pre.'jented  as  required  in  each  subregion. 

The  regions  chosen  for  this  specific  purpose  all  have  topography  which  is 
relativelv  well  represented  at  JO  km  resolution.  In  each  ca.se.  regions  are  defined  bv  the 

i'.  i  .  I  h.  ? .  ;Mf‘ 
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presented  by  showing  the  observ  ed  precipitation  across  a  topographic  cross-section  which 
is  chosen  based  on  observational  density  as  well  as  topographic  gradient.  Observations 
ate  then  compared  to  the  u  qf  derived  quantities.  Topography  requiring  finer  resolution  to 
be  represented  does  not  fare  well,  as  smoothing  tends  to  destroy  the  forcing  maxima  (a 
good  example  is  the  Santa  Cruz  mountains).  As  much  as  possible,  limitations,  as  well 
as  strengths,  corresponding  to  specific  regional  correlation  will  be  identified. 

2.  Topographic  Subregions 

5on  Pernar/finolSoti  Gdbricl 

This  subregion  was  particularly  hard  hit  by  storms  early  in  the  period. 
Figure  20  shows  the  regional  topography  to  consist  of  a  broad  flat  area  bordered  to  the 
northeast  by  the  abrupt  rise  of  the  San  Bernardino  range.  Figure  21  shows  the  actual 
obseived  precipitation  accumulations  along  with  the  derived  quantities,  for  0600  UTC 
-1800  IITC  01  March  1991  combined  with  a  representative  cross-section  of  the  regional 
topography.  This  particular  cross-section  is  representative  of  the  general  correlation 
within  this  subregion,  and  was  chosen  since  it  represents  the  period  of  maximum  rainfall. 

Tbe  wqf  parameter  seemed  to  do  a  particularly  good  job  heindling  the 
Mojave  Desert  stations  throughout  the  period,  while  showing  only  moderate  success  in 
the  upslope  regions.  This  may  best  be  explained  by  the  coarse  resolution  of  the 
topography  southwest  of  the  crest  Our  technique  sees  only  a  broad,  continuously  sloping 
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Figure  19.  Wqf  vs.  observed  precipitation  (inches,  vertical  axis)  for  San 
Bemardino/San  Gabriel  subregion  (0000  UTC  28  February  1991,  top)  and  Southern 
Sierra  Nevada  subregion  (1200  (.TTC  01  March  1991.  bottom) 
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region,  whereas  many  reporting  stations  lie  in  dissected  valleys  (i.e.  Riverside  at  301  m 
is  assigned  an  elevation  of  550  m  by  the  40  km  topography).  Also  note  (stations  1,  2  and 
3  on  Figure  21)  the  inflated  ■^  qf  values  close  to  the  ocean.  This  results  from  the 
smoothing  process  described  earlie*  . 

With  some  notable  exceptions  (i.e.,  Lytle  Creek  Ranger  Station  (LYCR 
#11 )).  does  a  decent  job  of  representing  the  regional  precipitation  showing  a  10% 
increase  in  the  correlation  coefficient  over  the  state-wide  case.  Forcing  responsible  for 
the  anomalous  observ^ations  (i.e.,  LYCR  and  CRES)  is  certainly  of  too  small  a  scale  to 
be  accurately  accounted  for  by  this  technique.  Still,  the  ability  of  wqf  to  correctly 
produce  a  specific  area  of  heavy  precipitation  (>1.5"  in  12  hours),  as  weU  as  define  the 
lesser  amounts  observed  in  the  Mojave  Desert  (stations  VICT,  BORO,  DAGG  and  MOJA) 
is  encouraging.  By  comparison,  NGM  12  hour  precipitation  totals  forecast  for  this  region 
were  all  less  than  0.4  inches  with  no  upslope  or  lee  structure  in  the  NGM  precipitation 
(not  shown). 

Of  course,  wqf  has  the  benefit  of  being  statistically  fitted  to  actual 
observations,  vvhere  the  NGM  forecast  represents  a  strictly  physical  quantity. 
Nonetheless,  the  statistical  fit  achieved  by  wqf  is  derived  objectively,  and  is  intended  to 
be  applied  in  all  conditions.  Srdl,  an  application  of  this  method  to  independent  data  is 
needed  to  determine  predictive  skill. 
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Fig:w»‘€  20-  Topography  (meters)  and  station  names  for  San  Bemardino/San 
Gabriel  subregion.  Precipitation  amounts  (inches)  axe  for  12  hours  ending  1800  UTC 
01  March  1991. 
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Figure  21.  Topography  (solid  line,  meters),  wq/' precipitation  (dashed,  inches), 
and  observed  precipitation  (points,  inches)  for  San  Bemardino/San  Gabriel  subregion 
(12  hour  period  ending  1800  UTC  01  March  1991). 


It  should  be  mentioned  that  the  significant  rain  events  in  both  this 
subregion  and  the  Santa  Barbara  subregion  which  follows  ended  by  1800  UTC  02  March 
1991.  As  a  result,  the  conditions  which  lead  to  deterioration  of  correlation  values  late  in 
the  period  were  not  observed  in  these  regions. 

b.  Santa  Barbara 

The  Santa  Ba.rbara  subregion  (Figure  22>  is  characterized  by  particularly 
good  ii'c^rf/precipitafi”'’  '.'M'.  l.Tii<'ii  'Oirni.  >  . .  -n  -.'  .  -.i  .  i.-ingiiir  hcni 
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-0  52  to  -0.85  for  the  dif^'erent  analysis  times.  Note  that  correlation  during  the  post¬ 
frontal  convective  period  of  1200  LfTC  28  February  1991  -  0000  UTC  01  March  1991 
shows  less  deterioration  in  this  region  than  in  the  state-wide  case  (-0.52  vs.  -0.22). 
Similar  to  the  previous  subregion,  lee  stations  (NE)  are  represented  highly  accurately, 
with  the  majority  of  the  variance  seen  in  the  opslope  regions  (SW).  Again,  subgrid-scale 
topography  likely  produces  large  variations  over  small  distances. 
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Figure  22.  Topography  (meters)  and  station  names  for  Santa  Barbara  subregion. 
Precipitation  amounts  (inches)  are  for  12  hours  ending  1800  UTC  01  March  1991. 


A  good  example  of  an  error  caused  by  insufficient  resolution  of  the 
topography  is  evident  at  time  0000  liTC  28  February  1991-  Figure  24  shows  the  wqf 
maximum  to  be  "transposed"  to  the  southeast  from  the  observed  maximum,  therefore 
limiting  co»Te)c>f'^''’ l'”  ''V- .  '  r  -;:' 'pb'  mthi  ; 
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Figure  23.  Topography  (so)id  line,  meters),  precipitation  (dashed,  inches), 
and  observed  precipitation  (points,  inches)  for  Santa  Barbara  subregion  (12  hour 
period  ending  1800  LUC  01  March  1991. 


region  (not  shown)  depicts  a  narrow  northeast  trending  ridge  in  the  vicinity  of  Thousand 
Oaks  (station  LECH),  northwest  of  which  is  the  Ventura  Valley  (station  SIMI).  This 
topography  is  represented  as  a  continuous  slope  at  40  km  resolution,  however.  As  such, 
the  strongest  orographic  ally  forced  component  of  wqf  is  located  over  the  Ventura  Valley 
(near  SIMI),  whereas  in  reality  this  forcing  maximum  is  located  between  the  Ventura 
Valley  and  Pine  Mountain  (station  PINM).  As  the  southeasterly  winds  become  more 
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Fig\ir«  24.  Wqf  and  observing  stations  with  12  hour  precipitation  totals  ending 
0600  UTC  28  February  1991  (inches)  for  the  Santa  Barbara  subregion  . 

favorably,  as  the  large-scale  winds  interact  with  topography  that  is  better  represented  at 

the  40  km  resolution. 

Overesttmation  of  precipitation  at  the  coast  (Figure  23.  stations  1-5)  does 

not  Sf'F'rn  Hp  fll*  ■  .  ■’  ’  ■  ''  -f  llv 
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field  still  prorluces  upward  orographic  forcing  well  into  oceanic  regions,  the  rapid  increase 
in  elevation  onshore  produces  orographically  enhanced  observed  precipitation  amounts 
eveti  near  the  ocean.  Moreover,  the  apparent  high  coastal  correlation  is  more  likely  a 
coincidence  as  opposed  to  an  increa.se  in  skill  over,  for  example,  the  San  Bernardino 
region.  As  with  the  previous  sub-region,  no  significant  precipitation  was  observed  after 
1800  Lire  02  March  1991. 

C-  Southern  Sierra  S'^rada 

The  most  notable  feature  seen  ijr  this  region  (Figure  25),  commonly 
kn<.>v,  Ti  as  the  High  Sierra,  is  the  lack  of  observations.  Still  we  can  make  at  least  one 

important  observation; 

*  VlVjrf  generally  overestimates  precipitation  for  these  inland  locations. 

This  trend  is  evident  tlrroughout  the  study  period  and  is  represented  m 
Figure  26.  We  can  generally  note  that  responsible  for  a  specific  rain  amount  along 
the  coast  (i.e.,  Santa  Barbara)  represents  a  somew'hat  lesser  amount  at  more  inland 
locations  This  feature  is  almost  certainly  a  function  of  the  available  moisture,  such  that 
inland  locations  require  a  greater  yyqf  to  maintain  saturation  than  is  required  near  the 
const  No  effort  is  made  to  quantify  this  relationship,  however  it  must  be  considered 
when  usijrg  ^yqf  based  on  inter-regional  statistics  (similai  to  adjusting  Model  Output 
Statistic  (MOS)  precipitation  forecasts  for  a  specific  location t 


Figure  25.  Topography  (meters)  and  station  names  for  Southern  Sierra  Nevada 
subregion  with  12  hour  precipitation  totals  (inches)  ending  1800  UTC  01  March 

1901. 

Note  that  Figure  2f'<  shows  and  observed  precipitation  both  before  and 
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Figure  26-  Same  as  in  Figure  23  except  for  Southern  Sierra  Nevada  subregion 
(12  hour  period  endmg  1800  UTC  01  March  1991.  top.  and  0600  I.TTC  04  March 
l’^91.  bottom) 
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orographic  effects  are  more  evklenr,  especially  considering  the  rough  equivalency  of 
precipitation  totals  at  the  San  Joaquin  Valley  station  Corcoran  (CORC)  and  amounts 
recorded  higher  in  the  Sierra.  The  low  slope  value  (-0.004)  of  the  linear  regression 
supports  this  assignment. 

rf.  Northern  Sierra  Nevada 

This  subregion  (Figure  27)  depicts  the  gradually  sloping  region  of  the 
Siena  Nevada  between  Sacramento  and  the  Lake  Tahoe  basin.  The  chosen  cross-section 
(Figure  2S)  coincides  roughly  w'ith  the  American  River  drainage,  which  is  represented 
particularly  well  by  the  ktn  topographic  resolution.  As  a  result,  we  see  very  good 
M  (7//precipitation  conelation  in  this  sub-region  prior  to  1200  LTC  02  March  1991,  with 
correlation  coefficients  between  -0.70  and  -0.80.  As  before,  we  note  the  overestimation 
of  piecipitation  in  the  Sacramento  valley  (WHEA),  as  well  as  a  general  tendency  to 
underestunate  precipitation  (slope  values  average  -0.006  in  this  subregion).  The  generally 
good  ftt  seen  in  the  first  period  at  the  summit  and  lee  stations  of  Soda  Springs  (SODA) 
and  Truckee  (TRUC)  respectively  represent  the  ability  of  this  technique  to  account  for  lee 
"spillover",  due  to  the  statistical  correlation  producing  a  positive  y-intercept. 

Perhaps  the  best  example  of  subgrid-scale  processes  affecting  the 
distribution  of  precipitation  is  seen  ,-jt  time  1200  UTC  04  March  1991  (Figure  27).  Note 
the  extraordinarily  high  precipitation  observations  over  Dormer  Summit,  including 
Tmckee,  as  compared  to  precipitation  estimated  by  nv/f  (regional  correlation  coefficient 
a»1'  I  sl'^P'"  \  nhiC"  ni  «■  ''  I  ' '  I  n  MM  ■  , ,  I .  . .  ■  . ' .  mmI  \  ;  I  :  i  iV)t 
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exoiiip'^  into  qviestiop  ny//'?  level  of  usefulness  given  this  type  of  svnoptic 

situation 

Of  particular  note  concerning  this  region  is  the  ability  to  compare  results 
obtained  here  to  those  obtained  over  a  similar  region  by  Colton's  (1976)  4.3  km  model. 
Though  Colton's  (1976)  use  of  an  hourly  precipitation  rate  prohibits  direct  comparison 
to  tlie  1 2  hour  rate  presented  here,  v-e  can  make  some  qualitative  assessments.  Other  than 
Colton  showing  considerably  greater  detail  within  the  upslope  region,  including 

the  sec<'tidaiy  ma;iimutn  produced  by  Mt.  Rose  rising  to  the  east  of  Lake  Tahoe,  wqf 
produces  a  similar  enhancement  pattern,  prior  to  1200  UTC  02  March  1991,  in  terms  of 
the  location  of  ma;timum  as  well  as  tlie  intensity.  Of  further  note  is  that  even  at  the  4.3 
km  grid  scale.  Colton's  (1^76)  calculations  at  lee  stations  generally  underestimate 
precipitation,  even  with  a  ''spillover  "  temi  included.  However,  Colton's  (1976)  fine  mesh 
resolution  allows  a  better  representation  of  convective  forcing.  In  these  types  of 
simations,  is  fat  inferior  to  Colton's  (1976)  model. 
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Figure  27.  Topography  (meters)  and  station  names  for  Northern  Sierra  Nevada 
subregion  with  12  hour  precipitation  totals  (inches)  ending  1800  UTC  04  March 
1091, 


The  final  sub-region  (Figure  29)  confains  not  oplv  the  most  complicated 
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FigvTE  29-  Same  ac'  in  Figvire  2?  except  for  Northern  Sierra  Nevada  subregion 
( 12  hour  accumulation?  ending  0600  blTC  02  March  1991  ,  top,  and  1800  UTC  04 
March  1991.  bottom). 
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02  March  1991  with  corresponding  correlation  coefficients  of  -0.68  and  -0.46  respectively. 
Also,  slope  values  for  this  subregion  average  -0  012.  which,  being  relatively  high,  is 
consistent  with  its  more  coastal  location.  Despite  the  obvious  smoothing  of  topographic 
features  at  this  resolution,  the  overall  wqf  representation  of  precipitation  is  quite 
encouraging,  in  particular  the  secondary  maximum  seen  between  Happy  Camp  (HAPP) 
and  Elk  Valley  (ELKV ).  Even  more  interesting  is  the  less  obvious  loss  of  correlation  late 
in  the  period.  This  may  be  due  to  this  subregion’s  northernmost  location.  Here, 
svnoptically  derived  subsidence  is  less  than  that  at  more  southern  locations  at  this  time, 
allowing  to  "create"  upward  motion  by  offsetting  the  synoptic-scale  subsidence. 
Clearly,  the  improved  correlation  is  only  coincidental  since  convective  processes  are  not 
accounted  for  in  any  region.  Another  possibility  is  more  consistent  forcing.  The 
pro.ximity  of  this  subregion  to  the  surface  low  pressure  systems  late  in  the  period  may 
have  resulted  in  a  more  consistent  low-level  southerly  flow.  In  this  situation,  though 
convection  is  still  unaccounted  for.  the  synoptic-  scale  winds  are  probably  acceptable 
representations  of  the  low-level  flow  in  this  subregion.  This,  combined  with  the  strong 
tendency  for  the  north  end  of  the  Sacramento  Valley  to  channel  winds  with  a  southerly 
component  may  have  allowed  urjrf  to  reach  a  level  at  least  of  the  same  order  of  magnitude 
a":  the  convective  forcing.  As  such,  this  period’s  correlation  is  likely  more  an  effect  of 
the  particular  local  situation  interacting  with  the  synoptic  scale,  rather  than  being 
more  accurate  in  this  subregion 
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Figure  29.  Topography  (meters )  and  station  names  for  Shasta-Siskyou  subregion, 
with  12  hour  precipitation  totals  (inches)  ending  1800  UTC  01  March  1991. 
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f.  Non-monntoinoti^  Kc^ions 


Though  this  study  is  aimed  at  representing  mesoscale  precipitation  within 
mountainous  regions,  it  is  also  hoped  that  including  the  effects  of  local  topography  will 
improve  distribution  in  low  lying  regions,  or  regions  in  which  topography  is  sub-grid 
scale  such  as  the  Santa  Cmz  Mountains  (Figure  31).  Without  going  into  great  detail,  we 
can  qualitatively  assess  the  success  of  vt-^f  as  weaker  where  topography  is  not  as  well 
defined.  As  before,  regions  of  less  dramatic  topography  are  most  susceptible  to  sub¬ 
synoptic  processes  dominating  the  local  distribution  of  precipitation.  This  is  especially 
evident  when  surface  forcing  is  not  represented  well  by  the  synoptic  analysis  (i.e.  upper 
level  features  and  convection). 

As  fa.r  as  mountainous  regions  which  are  poorly  represented  at  40  km 
resolution,  there  is  sufficient  evidence  to  suggest  improving  topographic  resolution  in 
order  to  identify  wqf  forcing  in  these  regions.  Again  we  refer  to  the  precipitation 
distribution  within  the  Santa  Cruz  mountains  (Figure  32).  Present  resolution  does  not  see 
the  steep,  narrow  nature  of  this  range,  but  observations  support  strong  local  orographic 
influence.  .Increasing  topographic  resolution  may  provide  sufficiently  large  quantities  of 
uv/f  to  "survive"  the  smoothing  process  realizing  that  this  representation  wtU  only  be  of 
use  if  the  synoptic  analysis  is  at  least  close  to  reality,  which  is  certainly  far  from 


guaranteed. 


Figure  31.  Topography  (meters)  and  station  names  in  the  vicinity  of  the  Santa 
Cruz  Mountains  with  12  hour  precipitation  totals  (inches)  ending  1800  UTC  28 
February  1991 

3.  Compari5on  to  NGM  Precipitation  Forecasts 

Earlier  we  established  the  rough  equivalency  of  the  NGM  precipitation  forecast 
with  the  term  used  here  to  represent  svnoptic-scale  forcing.  u  Of.  We  have  seen  examples 
of  the  improvement  the  topogiaphic  forcing  term  has  on  the  uOf,  field  in  representing 
ohser^ed  precipitation  during  the  first  period.  This  is  essentialh  the  same  level  of 
improve  me 'U  uvyi'  iH  i  u  .  i ,  ii,.--  i  i'..  -n  •  >  i 


Figure  32.  Wqf  9n6  observing  stations  with  12  hour  precipitation  totals  ending 
1800  28  February  1991  (inches)  in  vicinity  of  Santa  Cruz  Mountains. 


into  wqf).  Note  that  while  the  NGM  and  M-Ofarc  able  to  define  general  precipitation  areas 
reasonably  well,  there  is  virtually  no  structure  other  than  that  associated  with  the  synoptic- 
scale  fojcing  mazvima  (usually  coincident  with  700  mb  vertical  velocity),  and  therefore 
weak  correlation  w'ith  observations 

Late  in  the  period.  NGM  represents  precipitation  v.  ith  only  slightly  less  overall 
accuracy  than  yyqf.  Recall  that  these  times,  the  NGM  vertical  velocity  forecast 
predicted  upw-ard  motion  in  areas  which  verified  to  be  subsident.  Still,  its  convective 

pa.raine^^r  m;n  bp  . . miff  ( . in.  ..'  '  .  ...  ... 


VI  CONCLUSIONS 


A.  ArpLICATION  TO  FORECASTING 

It  i?  well  known  that  precipitation  amounts  can  var>'  dramatically  over  a  small  area. 
Inclusion  of  topography  into  fine-mesh  models  has  shown  favorable  representation  of  this 
variability  (Colton  1976.  Collier  1976).  Unfortunately,  the  data  network  needed  to  mn 
these  iiK'dels  is  not  generally  available  to  the  forecaster,  who  can  orrly  combine  local 
knowledge  vvitb  the  svnnptic-scale  piodncts  he  receives.  In  this  thesis  we  have  attempted 
to  provide  a  quantitative  estimate  of  mesoscale  precipitation  distributions  over 
mountainous  terrain,  without  requiring  mesoscale  information,  other  than  topography. 
Conclusions  can  be  summarized  by  the  following  points: 

-  C'iven  that  synoptically  mduced  forcing  is  well  represented  in  the  analyzed  or 
forecast  fields,  can  predict  a  substantial  degree  of  the  mesoscale  precipitation 
distribution. 

■  Most  of  the  discrepancies  between  wqf  and  the  observations  occur  within  upslope 
regions  (as  opposed  to  lee  stations)  due  to  unresolved  small-scale  topographic 
forcing. 

-  Correlation  is  greatly  unproved  when  considering  well  defined  topographic  regions. 

■  shows  significant  improvement  over  the  NGM  precipitation  forecast  fields. 

Coincident  with  our  conclusion  that  wqf  is  a  good  indicator  of  precipitation  under 

well  dr‘(iu''H  ,':vn''’pti'  -  .'I-  . . .  '.  - -  ,  ..i  I, 
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Tlterefo'e.  to  be  of  use  to  the  forecaster,  iie/she  must  first  determine  the  applicability  of 
u  (/f  based  on  the  synoptic  situation.  There  is  sufficient  evidence  to  suggest  wqf  may  be 
of  some  value  even  when  no  synoptic -scale  forcing  is  indicated.  In  this  case,  the 
forecaster  must  weight  the  field  accordingly,  but  be  aware  that  areas  of  consistently 
strong  Cow  may  still  be  represented  by  M  precipitation  maxima.  As  such,  assuming  wqf 
is  produced  as  a  forecast  field,  the  forecaster  must  identify  its  credibility  based  on  the 
following  parameters: 

-  Fresence  of  a  well-defined  thickness  ridge  ahead  of  a  front. 

-  The  potential  for  unstable  precipitation. 

It  these  conditions  can  be  correctly  evaluated,  wqf  will  be  sufficiently  useful  to  the 
forecaster  in  predicting  locally  heavy  precipitation  to  warrant  its  calculation  by  regional 
forecast  offices 

R.  RECOMMENDATIONS 

The  most  immediate  recommendation  would  be  to  incorporate  actual  NGM  forecast 
fields  into  the  wqf  calculation.  Tliis  would  provide  a  basis  for  verifying  the  wqf 
precipitation  against  the  actual  NGM  forecast  precipitation  amounts.  Such  an  application 
will  be  able  to  more  definitively  evaluate  wqf  ss  a  forecast  parameter. 

Secondly,  there  may  be  some  benefit  to  computing  over  even  finer  topographic 

resr'luti'''n  li  e  .  2'*  ■  i'’  if>>  tn- .  ii.  i  .  •  •..  i,  .i.  ..  l  in  he 
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improve  wqf  representation  of  precipitation  by  pinpointing  areas  of  maximum 
etrhancement  and  by  reflecting  the  orographic  influence  of  steep  topography  with  less 
horizontal  extent.  Tliere  will,  however,  most  certainly  be  a  point  at  which  no  increase  in 
wqf  accuracy  is  seen  with  increasing  topographic  resolution.  Furthermore,  objective 
restrictions  potentially  can  be  applied  to  the  calculation  of  wqf  \n  flat  regions  which  are 
adjacent  to  significant  topographic  barriers.  The  intent  here  is  to  limit  the  forcing  to  the 
actual  topographic  slope  ratlrer  than  distributing  it  over  a  larger  region  as  done  by  the 
ciinent  smoothing  process. 

A"  a  final  note.  1  wish  to  etnphasize  that  no  new  ground  has  been  broken  here  as 
far  as  o'u  understanding  of  the  precipitation  process.  I  have  simply  adapted  seeder-feeder 
theory  of  orographicalb  enhanced  precipitation  to  a  fomiat  which  only  requires  synoptic 
scale  infonnation  Thougli  the  results  are  far  from  ideal.  I  think  the  technique  exhibited 
sufficient  success  to  wanani  its  consideration  by  the  forecaster. 
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